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We present for the first-time the photonuclear cross section of calculated within the ab-initio No 
Core Shell Model framework. Realistic two-nucleon (NN) chiral forces up to next-to-next-to-next- 
order (N3LO), which have been softened by the similarity renormalization group method (SRG) to 
A = 2.02 fm“^, were utilized. The electric-dipole response function is calculated using the Lanczos 
method. The effects of the continuum were accounted for by including neutron escape widths derived 
from R-matrix theory. The calculated cross section agrees well with experimental data in terms of 
structure as well as in absolute peak height, (Tmax = 4.85 mb at photon energy oj = 23.61 MeV, and 
integrated cross section 85.36 MeV-mb. We test the Brink hypothesis by calculating the electric- 
dipole response for the first hve positive-parity states in ^°B and verify that dipole excitations built 
upon the ground- and excited states have similar characteristics. 


Electric-dipole transitions are an important excitation 
mode characterizing many facets of nuclear structure. 
Of particular interest is their strongly collective nature, 
which is manifested in what is known as the giant-dipole 
resonance (GDR). The GDR is ubiquitous in nuclei, and 
since its initial observation [1], much experimental and 
theoretical effort has been devoted to understanding its 
properties. The centroid of the GDR generally scales as 
the inverse of the nuclear radius, and experimentally is 
found to be ^ MeV, while the width is of the or¬ 

der 5 MeV [2]. Early on, phenomenological models were 
proposed by Goldhaber and Teller [3] and Steinwedel and 
Jensen [4] based on proton-neutron fluids that were able 
to describe the energy of the resonance. The width, on 
the other hand, was postulated to be due to the GDR 
damping into other nuclear modes of motion [5, 6]. A fur¬ 
ther, intriguing property is that a collective dipole mode 
exists on each state of the nuclear system, as hypothe¬ 
sized by Brink [7]. The resonant part of the photonuclear 
cross section has been calculated with semi-realistic in¬ 
teractions for "^He [8, 9], ®He and ®Li [10, 11] and ^Li 
[12]. More recent calculations for "^He utilizing modern 
realistic two- and three-body interactions have also been 
performed [13, 14]. These calculations have either used 
the framework of the hyperspherical harmonics (HH) ex¬ 
pansion [15, 16] or the No Core Shell Model (NCSM) [17- 
19]. Recently the Lorentz integral method [20, 21] was 
used in conjunction with coupled-cluster calculations to 
calculate the giant dipole resonance [22]. 

In this letter, we report on a theoretical study of the 
properties of the GDR for within the framework of 
the ab initio No Core Shell Model (NCSM) [17-19]. We 
study the convergence properties of the GDR as a func¬ 
tion of the model space size and present the photo-nuclear 
absorption cross section for ^°B. We demonstrate the in¬ 
fluence of more complex modes on the damping of the 
GDR as well as the influence of the neutron escape width 
on the dipole response. Finally, we test the Brink hypoth¬ 
esis by calculating the dipole response on positive parity 


excited states in ^°B and find a robust GDR built on each 
of these states exhibiting remarkably similar properties. 

The NCSM is a bound-state technique appropriate 
for light nuclei that uses as input realistic two- and 
three-body nuclear interactions. The NCSM determines 
the eigenenergies and wave functions of the nucleus by 
expressing the translationally invariant Hamiltonian in 
terms of antisymmetric combinations of single-particle 
harmonic oscillator (HO) states of frequency D. The size 
of the Slater determinant basis is determined by the to¬ 
tal HO quanta, A^maxj available in the system above the 
lowest configuration. Realistic interactions that make 
the link between structure and quantum chromodynam¬ 
ics explicit are derived using the effective field theory 
(EFT) for nuclear forces [23, 24]. In this work, we include 
nucleon-nucleon (NN) terms up to next-to-next-to-next- 
leading order (N3LO) [25]. To enhance convergence, an 
effective interaction was employed using the similarity 
renormalization group procedure (SRG) [26-29] with a 
momentum-decoupling value of A = 2.02 fm“^. At this 
A value, the binding energies of p-shell nuclei are repro¬ 
duced as though the calculation was performed with both 
the N3LO NN and N2LO NNN interactions [30-32]. In 
order to isolate effects of the strong interaction, we use an 
isospin-symmetric (isoscalar) interaction, and ignore the 
Goulomb interaction. All calculations were performed 
with Ml = 20 MeV. 

The dipole response function S{uj) on an initial state 
with angular momentum J and energy E is given by 
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where we assume the photon polarization is in the z- 
direction, the sum is taken over all initial orientations M 
and final states /, and w is the photon energy. is the 
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translationally-invariant dipole operator 
A 

= (1/2 - ZjA + i,,) , (2) 

i=l 

with the and fli denoting the proton number, single¬ 
particle radial and angular coordinates, respectively, and 
tz^i is the z-component of the isospin operator whose 
eigenvalues are ±1/2 when applied to a proton/neutron, 
respectively [33]. Finally, B{E1]J -A Jj) is the reduced 
transition probability. 

The response function, S{ui), was computed using the 
Lanczos algorithm as described in Refs. [34-39]. First, 
Dz is applied to an eigenstate and angular momentum 
J/ is then projected on this vector. This new vector is 
then used as an initial pivot for the Lanczos procedure, 
and n Lanczos iterations are performed giving n eigenval¬ 
ues with energy Ef. The reduced transition probability is 
then extracted from the overlap between the initial pivot 
and the Lanczos eigenvectors, and Eq. (1) is evaluated. 
We note that of these n Lanczos eigenvalues, only a sub¬ 
set of them is fully converged (generally just the extreme 
values), but with Eq. (1), these n Lanczos eigenvalues 
reconstruct the first 2n — 1 moments of S'(w). S{uj) was 
computed with n up to 300. All of the many-body calcu¬ 
lations, including the Lanczos calculation of the strength 
function distribution, were carried out using the BIG- 
STICK configuration-interaction code [40]. 

The El response function for the ground-state of 
is shown in Eig. 1 as a function of increasing basis size 
Amax- Overall, the GDR excitation spectrum shifts to 
lower energy with increasing Amax, reflecting a some¬ 
what slower convergence for these negative-parity excited 
states relative to the ^°B ground state. 



FIG. 1. (color online) The components of the strength func¬ 
tion for the ground-state {J^ =3^) as a function of in¬ 
creasing basis size Amax at = 20 MeV. The Amax = 5 — 9 
strength functions were calculated with 300 iterations. 


Fig. I illustrates one of the principal broadening mech¬ 


anisms for the GDR. In particular, as Amax increases, the 
simple collective 1-particle, 1-hole (Ip-lh) nature of the 
GDR couples to more complex states, and the strength 
“spreads out” across a wider range of states. In partic¬ 
ular, the GDR is “damped” via mixing with these more 
complex low-lying 3p-3h (and higher) states. 

A second broadening mechanism must be addressed 
in order to compare with experimental data. This arises 
from the fact that most of the states excited by the dipole 
operator lie above the threshold for particle emission; in 
particular neutrons at Sn = 8.436 MeV [41]. Thus, the 
dominant decay mode of these states is neutron emission, 
which occurs on a short time scale. We can estimate 
the width of these excited states with the neutron-escape 
width, r„, from the R-matrix theory 

r„(F;,) = 272p^p,0/, (3) 

I 

where E^ is the excitation energy, = {hc )"^/is 
the single-particle width, fi = ((A — 1)/A)tojv is the re¬ 
duced mass of the A-nucleon system in terms of the nu¬ 
cleon mass TOjv, and the radius is i? = 1-2[(A — 1)3 ± 1]. 
The penetration factor Pi = knRVi depends on the wave- 
vector kn = \/2^En {En is the neutron kinetic energy) 
and the penetrability, Vj, of the outgoing neutron with 
orbital angular momentum 1. Since ^°B is a p-shell nu¬ 
cleus and the GDR is dominated by Ip-lh excitations 
from the Op shell into the Od — Is shell, neutron emis¬ 
sion of these states will be dominated I = 0 and I = 2 
partial waves. For s-waves, Vq = 1, while for d-waves, 
V2 = (A:i?)V(9 ± 3(kPy + (kR)‘^) [42]. 

The relative strength for each partial wave in Eq. (3) 
is given by the spectroscopic factor, 

01 = [(9R; JdTdl X (n; {nlj-, 1/2)|] \^^B- JfTf), (4) 

which is a measure of the overlap of the neutron±®B 
with the neutron-unbound ^°B (the neutron is in the HO 
orbital state {nlj) with isospin t = 1/2). The magni¬ 
tude of the width is strongly governed by the outgoing 
kinetic energy (^ y/En), which in turn is governed by 
the distribution of the spectroscopic strength to states 
in ®B {En = Ex — Sn — Ed). Indeed, in the simplest 
approximation, where the entire spectroscopic strength 
is collected in the ground state of ®B, r„ increases dra¬ 
matically with excitation energy, and is unphysical. On 
the other hand, due to the large dimensions present it is 
also not practical to calculate the spectroscopic strength 
in the largest model spaces and for highly excited states. 
Instead, we build a model for the spectroscopic strength, 
Qf{En + Sn), based on fairly complete calculations for 
Ainax = 3 and 5. The spectroscopic strength is modeled 
with a normalized Gaussian that reproduces the total 
strength as well as the first and second moments of the 
calculated A„iax = 3 and 5 results, which is feasible up 
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to Ex ~ 30 MeV. The neutron-escape width is then esti¬ 
mated by 

T^{Ex) = V / " " PiQfiEn + S^)dEr,. (5) 

l -^0 

In Fig. 2, we compare the actual widths calculated by 
Eq. (3) with those of the model of spectroscopic factors 
from Eq. (5). We find that r„ increases with excitation 
energy, but much less dramatically than if the decay oc¬ 
curred to the ground state. Further, the width appears to 
saturate to a value of 10-12 MeV at « 40 — 50 MeV. 
This is largely due to the fact that the distribution of 
spectroscopic strength occurs to states with one less os¬ 
cillator quanta than the initial state, and thus tends to 
cap the outgoing kinetic energy at roughly Ml. For what 
follows, we compute the neutron escape width using the 
model for the spectroscopic strength constrained by the 
.Vnax = 5 calculation. 



FIG. 2. (color online) The widths of negative parity states 
as a function of the excitation energy {Ex) above the neu¬ 
tron separation energy {S„) are shown by the blue (solid) line 
when the spectroscopic factors are taken into account as in 
Eq. (3) whereas the red (dashed) line shows the result of us¬ 
ing a model of the spectroscopic factors as in Eq. (5). The 
calculation is performed at Vmax = 5. 


The photo-absorption cross section in the unretarded 
dipole-approximation is related to the strength function 
by (j{uj) = d-TT^aw {S{uj) — S{—oj)), where a = |^ is the 
fine-structure constant [2]. For states above the neu¬ 
tron separation energy, the discrete lines in S{uj) were 
replaced with a normalized Breit-Wigner function cen¬ 
tered at the eigen-energy Ef with width r„(E/). The 
photo-absorption cross section is then 


ct(w) = dTT^a SfEf 


'r„(Ey) 




(ccr„(£;/))2’ 


( 6 ) 


with Tn{Ef) given by Eq. (5) as shown in Fig. 2 and 
Sf = B{E\]J ^ Jf)/2>. The photo-absorption cross 


section calculated for is shown in Fig. 3 a), for 
-Vnax = 5 — 9 as well as an extrapolation to lYmax = 
oo . The vertical lines show the discrete spectrum at 
-Vnax = oo, while the continuous lines show the cross 
section with the neutron-escape width folded into the 
spectrum. The V^ax = oo result was obtained with an 
exponential extrapolation [43], where the absolute ener¬ 
gies of the ground state and the first large peak in the 
dipole-response function shown in Fig. 1 (see the peak 
at 20 MeV for Vmax = 9) for = 3 — 9 were fit to the 
function E{Nynax) = 4l exp(—BVmax) + C. The parame¬ 
ter C is then the energy at = oo. This procedure 

leads to a final excitation energy for the fitted dipole 
state of 18.209 MeV, which is 1.881 MeV lower than the 
-^max = 9 result. To estimate the = oo strength 

function, we assume this shift to be uniform across the 
entire spectrum, and shift the = 9 dipole response 
lower by 1.881 MeV. Overall, the cross section is seen to 
have a main peak at about 23 MeV for Vmax = oo, fol¬ 
lowed by a slightly smaller hump near 33 MeV. We also 
note that the main peak in the spectrum tends to shifts 
down in photon energy (w) as well as decreasing in height 
as the model space increases. 

In Fig. 3 b), we compare the V^ax = 9 and lYmax = oo 
calculation with experimental data. Hughes et al [44] 
and Ahsan et al [45] determine the cross section from 
bremsstrahlung sources, whereas Kneissl et al [46] used 
a quasi-monoenergetic photon beam from positron anni¬ 
hilation in flight. The bremsstrahlung data is in good 
agreement with each other whereas the Kneissl data is 
about 10% lower at the peak and has a slower fall-off in 
the tail of the cross section. In Ref. [46] it is argued that 
cross sections determined from bremsstrahlung sources 
have to be extracted from yield curves by an unfold¬ 
ing procedure, which may introduce spurious structure 
in the cross section. Furthermore, comparison between 
absolute cross sections for ®Li and ^Li measured at Liv¬ 
ermore [47, 48] using monoenergetic photons and those 
by Hayward [49] with bremsstrahlung sources are also 
about 10% lower. The peak in our calculated cross sec¬ 
tion is slightly higher in photon energy than experimen¬ 
tal data by about 1.5 MeV. The absolute peak height of 
CTmax = 4.85 mb at w = 23.61 MeV compares well with 
the Kneissl data and shows a similar pattern in the tail 
of the cross section. The N^ax = oo integrated cross sec¬ 
tion, a{ijj)du) = 85.36MeV-mb compares well with 

the experimental value of 83.1 ± 1.2 MeV-mb [46]. An¬ 
other evident feature is that both the experimental data 
and the theory calculation bear out a two-hump structure 
in the cross section implying that ^*^B may be deformed. 
At Vinax = 9, our calculated quadrupole moment for the 
= 3+ state is Q = 6.22fm^. Using a hydrodynam- 
ical model of the nucleus [50] one can estimate the de¬ 
formation parameter /32 ,ncsm ~ v^Q/( 3ZR2) ^ Q 37^ 
where we have taken our calculated radius as i? = 2.11 fm 
[51]. A similar calculation using the experimental data 
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of Q = 8.477 fm^ [41] and rms charge radius R = 2.58 fm 
[52] leads to /32,exp = 0.34. Varying the fit parameters of 
the spectroscopic factor model changes the peak-height 
of the cross section by about 10% but does not affect the 
threshold or tail region of the calculated cross section. 



FIG. 3. (color online) a) The photonuclear cross section as 
a function of V^ax = 5 — 9,oo. Note that the main peak 
slowly shifts down in photon energy (ta) as well as decreases 
in height. The energy-weighted strength, SfEf, is shown for 
Vmax = oo with Vertical lines, b) The calculated cross section 
at Nmax = C !0 (solid) is compared to three sets of experimen¬ 
tal photoneutron cross section data (symbols). Our absolute 
peak height of cTmax = 4.85 mb at ta = 23.61 MeV compares 
well with the Kneissl data. 


Lastly, we address the Brink hypothesis, which postu¬ 
lates that if the GDR is observed for the ground-state, 
then all excited states will similarly support a dipole res¬ 
onance [7]. In Fig. 4, we show results of the GDR photo¬ 
absorption cross calculated on the five lowest states in 
the NCSM spectrum for relative to the energy of the 
initial state. Our calculations give remarkably similar 
cross sections for each of the states, with a strong collec¬ 
tive mode at approximately the same excitation energy 
and width. We thus report the first fully microscopic cal¬ 
culation, based on realistic nucleon-nucleon interactions, 
that demonstrate the universal feature of the giant-dipole 
resonance in nuclei. 

In conclusion, we have presented the first NCSM cal¬ 
culations of the dipole response function of ^°B using 
realistic two-body chiral N3LO interactions that have 
been softened by SRG to A = 2.02 fm“^. To deter¬ 
mine the photonuclear cross section we assigned finite 
widths to the continuum-discretized spectrum by using 
the physics of neutron-escape widths from R-matrix the¬ 
ory. A model was built to mimic the behavior of the 
spectroscopic factors as a function of the excitation en¬ 
ergy of the ^°B particle-unbound states. The calcu¬ 
lated cross section agrees well with the experimental 
data in terms of structure as well as in absolute peak 



FIG. 4. (color online) The dipole response of the first five 
states in is shown. The dipole response of the excited 
states has been shifted down in photon energy by their re¬ 
spective excitation energy Ex,i- 


height CT„iax(23.61) = 4.85 mb and integrated cross sec¬ 
tion 85.36 MeV-mb. Our calculations support the Brink 
hypothesis and illustrate that excited states may have a 
different GDR shape than that of the ground-state. The 
methodology presented here represents the first step in 
creating a NCSM framework for calculating electroweak 
observables for light nuclei in general. 
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